Altered subgenomic RNA abundance provides unique insight into SARS-CoV-2 B.1.1.7/Alpha variant infections.

B.1.1.7 lineage SARS-CoV-2 is more transmissible, leads to greater clinical severity, and results in modest reductions in antibody neutralization. Subgenomic RNA (sgRNA) is produced by discontinuous transcription of the SARS-CoV-2 genome. Applying our tool (periscope) to ARTIC Network Oxford Nanopore Technologies genomic sequencing data from 4400 SARS-CoV-2 positive clinical samples, we show that normalised sgRNA is significantly increased in B.1.1.7 (alpha) infections (n = 879). This increase is seen over the previous dominant lineage in the UK, B.1.177 (n = 943), which is independent of genomic reads, E cycle threshold and days since symptom onset at sampling. A noncanonical sgRNA which could represent ORF9b is found in 98.4% of B.1.1.7 SARS-CoV-2 infections compared with only 13.8% of other lineages, with a 16-fold increase in median sgRNA abundance. We demonstrate that ORF9b protein levels are increased 6-fold in B.1.1.7 compared to a B lineage virus in vitro. We hypothesise that increased ORF9b in B.1.1.7 is a direct consequence of a triple nucleotide mutation in nucleocapsid (28280:GAT > CAT, D3L) creating a transcription regulatory-like sequence complementary to a region 3' of the genomic leader. These findings provide a unique insight into the biology of B.1.1.7 and support monitoring of sgRNA profiles to evaluate emerging potential variants of concern.

ACE2-Independent Interaction of SARS-CoV-2 Spike Protein with Human Epithelial Cells Is Inhibited by Unfractionated Heparin.

Coronaviruses such as SARS-CoV-2, which is responsible for COVID-19, depend on virus spike protein binding to host cell receptors to cause infection. The SARS-CoV-2 spike protein binds primarily to ACE2 on target cells and is then processed by membrane proteases, including TMPRSS2, leading to viral internalisation or fusion with the plasma membrane. It has been suggested, however, that receptors other than ACE2 may be involved in virus binding. We have investigated the interactions of recombinant versions of the spike protein with human epithelial cell lines that express low/very low levels of ACE2 and TMPRSS2 in a proxy assay for interaction with host cells. A tagged form of the spike protein containing the S1 and S2 regions bound in a temperature-dependent manner to all cell lines, whereas the S1 region alone and the receptor-binding domain (RBD) interacted only weakly. Spike protein associated with cells independently of ACE2 and TMPRSS2, while RBD required the presence of high levels of ACE2 for interaction. As the spike protein has previously been shown to bind heparin, a soluble glycosaminoglycan, we tested the effects of various heparins on ACE2-independent spike protein interaction with cells. Unfractionated heparin inhibited spike protein interaction with an IC50 value of <0.05 U/mL, whereas two low-molecular-weight heparins were less effective. A mutant form of the spike protein, lacking the arginine-rich putative furin cleavage site, interacted only weakly with cells and had a lower affinity for unfractionated and low-molecular-weight heparin than the wild-type spike protein. This suggests that the furin cleavage site might also be a heparin-binding site and potentially important for interactions with host cells. The glycosaminoglycans heparan sulphate and dermatan sulphate, but not chondroitin sulphate, also inhibited the binding of spike protein, indicating that it might bind to one or both of these glycosaminoglycans on the surface of target cells.

Subgenomic RNA identification in SARS-CoV-2 genomic sequencing data.

We have developed periscope, a tool for the detection and quantification of subgenomic RNA (sgRNA) in SARS-CoV-2 genomic sequence data. The translation of the SARS-CoV-2 RNA genome for most open reading frames (ORFs) occurs via RNA intermediates termed "subgenomic RNAs." sgRNAs are produced through discontinuous transcription, which relies on homology between transcription regulatory sequences (TRS-B) upstream of the ORF start codons and that of the TRS-L, which is located in the 5' UTR. TRS-L is immediately preceded by a leader sequence. This leader sequence is therefore found at the 5' end of all sgRNA. We applied periscope to 1155 SARS-CoV-2 genomes from Sheffield, United Kingdom, and validated our findings using orthogonal data sets and in vitro cell systems. By using a simple local alignment to detect reads that contain the leader sequence, we were able to identify and quantify reads arising from canonical and noncanonical sgRNA. We were able to detect all canonical sgRNAs at the expected abundances, with the exception of ORF10. A number of recurrent noncanonical sgRNAs are detected. We show that the results are reproducible using technical replicates and determine the optimum number of reads for sgRNA analysis. In VeroE6 ACE2+/- cell lines, periscope can detect the changes in the kinetics of sgRNA in orthogonal sequencing data sets. Finally, variants found in genomic RNA are transmitted to sgRNAs with high fidelity in most cases. This tool can be applied to all sequenced COVID-19 samples worldwide to provide comprehensive analysis of SARS-CoV-2 sgRNA.